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In a previous paper, a model for queuing processes with correlated 
inputs was analyzed. To illustrate the application of those results, we 
model a concentrator of slow terminals. The sources of data, the ter- 
minals, generate data at a slower rate than the output speed of the 
buffer. In special cases, we obtain closed-form expressions for the 
generating function of the equilibrium queue size distribution. For the 
general case, we describe a computational procedure to obtain the 
distribution and the average of queue size in steady state. The nu- 
merical results obtained using this procedure are presented for a family 
of problems in which each message consists of two packets separated 
by a fixed time interval. 

I. INTRODUCTION 

A data communications network may be constructed by connecting 
together terminals and switching nodes so that each terminal is con- 
nected to just one node and the nodes are connected together in a 
more-or-less redundant fashion. All connections are by means of 
transmission lines. Those between terminals and nodes are called access 
lines, while those between one node and another are called trunk lines 
(Fig. 1). For economy, the access lines commonly have smaller bandwidth 
than the trunk lines. 

The character of the traffic carried by a data network depends in part 
upon the type of terminal connected to it. Keyboard and display ter- 
minals transmit and receive data messages that are typically less than 
a few hundred characters in length. These terminals operate at speeds 
up to 1.2 Kb/s. Batch stations transmit and receive data in larger quan- 
tities and typically operate at speeds up to 9.6 Kb/s. Trunk lines, arid 
the connections to computers, typically operate at about 50 Kb/s. Thus, 
we find that traffic in a data network is not uniformly distributed either 
among the terminals or among the various transmission lines. 

Analysis of delay and the probability of queue overflow is most simply 
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SWITCHING NODE 
LOW-SPEED ACCESS LINE 
= HIGH-SPEED TRUNK LINE 

Fig. 1 — Network topology. 

obtained by assuming that all traffic is uniformly distributed between 
the terminals and that all transmission lines operate at the same speed. 
Using these assumptions, Chu 1,2 has studied two cases. In one case, equal 
size packets (individual characters perhaps) are generated randomly by 
the terminals. In another case, messages are generated at random, but 
each message consists of a random number of packets which all enter 
the network in one instant. The case of mixed input traffic was studied 
by Chu and Liang. 3 

Consider the situation when terminals randomly generate messages 
consisting of several fixed-size packets. The data are fed into a switching 
node over access lines that are substantially slower than the trunk lines 
used to carry data out of the node. Several packets may be transmitted 
on a trunk line in the time that it takes to transmit one packet on an 
access line (Fig. 2). Thus, packets of one message which are transmitted 
consecutively by a terminal will arrive periodically at the node, and the 
period will be greater than the period of packet transmissions on the 
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trunk line. Given that packet arrivals are not entirely random, one would 
expect the buffer storage requirements in the switch to be less than is 
indicated by Chu and Liang's analyses. 

In this paper, we study the behavior of a switching node that receives 
data from a (large) number of low-speed access lines. The data are re- 
ceived in the form of packets of a fixed size. As the packets arrive, they 
are placed in a buffer, which is a first-in-first-out queue. In an actual 
communications network, the buffer has a finite size, and a packet is lost 
if the buffer is full when it attempts to enter it. The buffer transmits 
packets at a uniform rate onto a high-speed trunk, provided the buffer 
is not empty. A crucial question is how large the buffer should be in order 
that the probability of packet loss should be less than 10~ 4 , say. In the 
mathematical analysis of the single queue corresponding to this buffer, 
we consider a buffer of unlimited size so that no overflow is possible, and 
we calculate the steady-state probability that the buffer content (i.e., 
the number of packets in the buffer, or queue size) exceeds the proposed 
size of the finite buffer. We refer to this quantity as the probability of 
overflow, since it is usually used to estimate the actual overflow proba- 
bility, when the probability of packet loss is very small. 

Some of the symbols used are listed in Section II. The mathematical 
model, which includes assumptions concerning the packet arrivals, is 
discussed in Section III. This model was analyzed by two of the authors, 4 
and formulas for calculating the equilibrium queue size distribution are 
summarized in Appendix A. These formulas involve some marginal 
distributions, and some polynomials which have to be determined. Ex- 
plicit analytical expressions for the coefficients in these polynomials are 
given for some particular examples in Section IV. The computation of 
the coefficients for more general examples is discussed in Section V. The 
computation of the marginal distributions is discussed in Section VI. 
Some numerical results are presented in Section VII. 

II. NOTATION 

b n buffer content at time n 

z n number of packets entering buffer in time interval 

(n,n + 1] 
k + 1 number of time intervals required for message arrival 

a) nonnegative integers with a' > 

(Xn, . . . ,x l n ) independent identically distributed vector of 

nonnegative integer valued random variables. 

III. MATHEMATICAL MODEL 

The behavior of a switching node is modeled by considering the state 
of the buffer at discrete times. Suppose that the buffer transmits one 
packet in a unit time interval. If b n denotes the buffer content at time 
n, the buffer content at time n + 1 is 
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b n +i = b n -l + z n ifb n >l 
= z n if b n = 



or equivalently 



>n+l 



= (b n -l) + + z n , (1) 



where z n is the number of packets entering the buffer in the time interval 
(n,n + 1], The peculiar feature of our queuing problem arises from the 
fact that the access lines are slower than the trunk line on which the 
output of the buffer is transmitted. For example, the 9.6 Kb/s access lines 
operate approximately at only one-fifth of the speed of a 50 Kb/s trunk. 
Then, in our discrete model, messages arriving from the access lines will 
consist of a random number of packets separated by five units of 
time. 

Let us first focus our attention on the case when there are two packets 
to each message. If the packets are separated by d units of time, then the 
number of packets entering the buffer in the interval (n,n + 1] will be 
equal to the number of first packets in messages arriving in that interval, 
plus the number of second packets in messages whose first packets ar- 
rived d intervals earlier. Then z n = x n + x n -d, where x n is the number 
of first packets arriving in the interval (n,n + 1]. Now consider the case 
in which the messages are transmitted by N independent sources. The 
probability is zero that the first packet of a message from a source enters 
the buffer in a unit time interval if the first packet of a message from that 
source entered the buffer in one of the previous (2d — 1) time intervals, 
and otherwise it is p. Then 

Pr{x n = i |*»-v = fr / = 1.2, • • • } = ( N ~ 7 )p I0 (l - p)"-'-'°, 

where / = Sj^J" 1 ij. Thus, x n and x n -d are not independent random 
variables. However, if N -*• °° with Np = X fixed, then Pr{x n = io) ""*" 

(e~ x \ io )/io\, independently ofx n -j, j = 1,2, Hence, if the number 

of sources is large, it is reasonable to assume that the random variables 
x n , n = 0,1, . . . , are independently and identically distributed, and we 
will not restrict ourselves to the Poisson distribution. 

Next we consider the case where each message consists of either two 
packets, separated by d units of time, or of only one packet. If we let x\ 
denote the number of first packets of two-packet messages arriving in 
the interval (n,n + 1], and x\ denote the number of single packet mes- 
sages arriving in the same interval, then z n = x\ + x\-d + jc„. In order 
to allow for randomness in the number of packets to a message (either 
one or two), x\ and x\ may be dependent on each other. For example, 
if the number of packets in a message is two with fixed probability 
1 — p, and one with probability p, where < p < 1, then E(tf nl t2 n2 ) = 
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6[(1 - p)ti + pt 2 ], where Q(t) is the distribution for the number of 
messages arriving in a unit time interval. As before, when there is a large 
number of independent sources, it is reasonable to assume that (x\,x%), 
n = 0,1, . . . , are independent identically distributed vector random 
variables. 

There are some obvious generalizations of the above special arrival 
processes. For instance, the random number of packets to a message may 
be as large as m ^ 2. Also, there may be slow access lines with several 
different speeds, so that the packets in a message may be separated by 
di units of time, i = 1, . . . ,r. This led two of the authors 4 to consider 
arrival processes of the form 

Zn = £ £ «}*U (2) 

i=l;'=0 

where the vector nonnegative integer valued random variables 
{(*£, . . . ,x l n )) are independent and identically distributed, and a) are 
nonnegative integers with a' > 0. In Appendix A we summarize the 
relevant results pertaining to the steady-state distribution of the queue 
size corresponding to (1) subject to (2). 

IV. EXPLICIT EXAMPLES 

To calculate the generating function of the steady-state queue size 
from (37) and (38), it is necessary to know the polynomials c r (s), r = 
1, . . . ,k, as well as the quantities <t> rv (s), r = 0,...,k. Hence, from (39), 
we need to know the constants c j. We now give explicit analytical results 
for some particular examples. 

We first consider the arrival process 

z n = x l n + x l n . k + x 2 n + xl (3) 

where k ^ 2, and 

E{t\h x Mh = e[(i - P )tx + P t 2 ]*(t3), (4) 

with < p < 1 fixed. This corresponds to arrivals from two different 
classes of sources. One class of sources sends messages which consist 
either of two packets separated by k units of time with probability 
1 - p, or of one packet with probability p. The other class of sources sends 
messages which consist of just one packet. From (2), (3), and (4), and the 
definition of v rn in (31), it follows that 

<t> ru (s) = 0(s)*(s), r=0,...,k-l, 

<t>kv(s) = 0[(1 - p)s 2 + ps]*(s). (5) 

It remains to give the values of c j, which we do for k = 2 and 3, with < 
p < 1, and for k = 4 with p = and ^(s) = 1. It is shown in Appendix C 
how to determine which constants c j occur in (46). It is also shown how 
the values of c j were calculated in the case k = 3. 
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Cm -— — -, ; ^j— — cqoo, (9) 



Let 

0(s) = £ p^, *(s) = £ qis 1 . (6) 

i'=0 i=0 

For example, the case p, = e~ x X l /i! corresponds to a Poisson distribution 
for the number of messages arriving from the first class of sources in a 
unit time interval. For k = 2 the nonzero coefficients c j are 

coo, en = (1 - p)pi<7oCoo, (7) 

with coo + Cn = 1. For k = 3 the nonzero coefficients Cj are 

cooo, Con = (1 - p)pi<7oCooo> 

C122 = (1 - P) 2 P?<70COOO, C 2 22 = U ~ p) 2 P0P2<?0C000, (8) 

and 

(1 - p)(? [pi(l +Pq(?i) + 2pp p 2 go] 

[1 - (1 - p)p Pi9o] 

with 2cj = 1. 

For fe =4we take p = and ^(s) s 1, corresponding to z n = x n + x n -^ 
with E(s Xn ) = 0(s). In this case, there are 14 nonzero coefficients Cj, 
namely 

coooo = Co, coon = Pico, C0122 = P1C0, 

C1233 = P?C0, C0222 = P0P2C0, C3333 = P0P3C0, 

C1333 - C2333 = C2233 = P0P1P2C0, (10) 

and 

com = PiAco, C1122 = PiAc , 

C1222 = Pi(l + PoPi)Ac , C2222 = P0P2U + PoPi)Ac , 

cim = Pi[l + Po(Pi + PoP2)(l + P0P1)] Ac , (11) 

where 

A = {1 - p pi[l + p 2 (pi 2 + PoP 2 )(l + PoPi)]}" 1 . (12) 

We now consider a class of arrival processes z n for which the polyno- 
mials c r (s), r = 1, . . . ,k, are, in fact, constants. Then the first two mo- 
ments of the equilibrium queue size distribution, Eyo and Eiyl), may 
be expressed, with the help of (59) to (61) and (63), in terms of the first 
three moments of u rn , since c' r (1) = and c" r (l) = 0. The class of arrival 
processes we consider corresponds to 

a) >0, j = 0, . . . J it i=l,... ,/, 

a) = 0, otherwise, (13) 
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in (2). We will show that y 0n = implies y rn = 0,r= I, . . . ,k, and hence 
that 

0(O,si, . . . ,s k ) = m = 1 - Ev kn . (14) 

This implies, from (39), that c r (s) = n, r - 1, . . . ,k. 

Now, since b n = y 0n , it follows from (29) and (13) that;yon = implies 
that 

*'„_,_! = 0, ; = 0, . . . J it i=l,...,l. (15) 

If r + 1 ^ ji, then a) ■ for j > r + 1. If, on the other hand, r + 1 < ;',-, 
then x|,_i+r - for r + 1 <;' < ji, and a) = for; > ;,. Hence, a)x l n - j+r 
= for; = r + 1, . . . ,k, i = 1, . . . ,/, r = 0, . . . ,k - 1. Thus, from (30),y n 
= implies y rn = 0, r = 1, . . . ,&, as was asserted. 

As a particular example, we consider messages which consist either 
of one packet or of two or more packets (up to a maximum number) 
which arrive in consecutive time intervals. Such is the case when the 
access lines have the same speed as the output line. For this example, 
the arrival process is 

k+i i-l 
2n=EI *U> (16) 

i=i y=o 
Hence, from (2), / = k + 1 and 

. = 1, ;' = 0, . . . ,i - 1, i = 1, . . . ,k + 1 (1?) 

; 0, ;' = i',... ,k,i = 1, . . . ,k 
so that (13) is satisfied, and c r (s) = n, r = 1, . . . ,k. From (31) it follows 
that 

/ij=min(i,r+ 1), r = 0,...,k, i = 1, . . . ,k + 1 (18) 

and 

k+i 

VOn = H x n, 
i=l 

r k+l 

V rn ='Lix i n+(r+l) L X l m T - 1, . . . ,k. (19) 

i=l i=r+l 

Let pi > 0, i = 1, . . . ,k + 1, be the probability that there are i packets 
in a message, where SfJj 1 p, = 1. Then, 

Ear ••• uu) = e(^z Piti)- (20) 

Hence, from (19), 

<f> 0u (s) = Q(s), 

<t>rv(s) = Q(£ Pis'+ E PfS^ 1 ), r-l,...,fc. (21) 
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V. COMPUTATION OF c r (s) 

As described in Appendix A [see (34)], the constants |cjJ determine 
the polynomials c r (s) from (39). In principle, using a program like AL- 
TRAN, 5 which performs symbolic manipulations, it is possible to sub- 
stitute (34) into (33) and equate coefficients of like powers on both sides 
of (33) to get the equations for the \cj\. However, we will describe an al- 
ternate method that was used for obtaining the numerical results pre- 
sented in Section VII. From (36) it is seen that 

f MC JV2 ■ ■ ■ jk ~ "ojij2 ■ ■ ■ ik- (22) 

The method presented here arrives at equations for {P j} directly from 
the equations for |Pi): from (32) 

" i = *- Pio~Ji.i\-J2, ■ ■ ■ ik—jk*i 
jo=0,jtA 

+ E Ph-h+U i-y'2+i ik-jh+iPy ( 23 ) 

jo>0,jtA 

The sequence of programs that were used to determine Pq/, . . j k are 
described briefly below. 

(i) The first step is to generate a sufficient number of equations from 
(23). Starting with i = (0,0, . . . ,0), the right-hand side of (23) is found 
in symbolic form. Using a test for determining admissible states, every 
index j appearing on the right-hand side of (23) that corresponds to states 
not communicating with (0,0, ... ,0) is omitted. Corresponding to each 
new index that arises on the right-hand side of (23), a new equation is 
generated from (23) by setting i equal to the new index. This process is 
terminated when every new index generated has io S> k. This whole 
process is repeated starting with each index with £q = 0, i.e., for each 
•Pq/i . . . jk- The total number of equations and the total number of un- 
knowns are counted. The output of this program is this set of equations 
in symbolic form. The number of equations turns out to be always less 
than the number of unknowns. However, we know 4 a set of linear ho- 
mogeneous equations must exist for P j 1 . . . j k . In the examples consid- 
ered, visual inspection revealed a few substitutions which made the 
number of equations one less than the number of unknowns. The details 
of a specific procedure for accomplishing this, in the case / = 1 in (2), have 
recently been given by Massey and Morrison. 6 

(ii) The normalizing constant which determines Po 71 . . j k uniquely 
is easily seen to be (m/SPq/i . . . j k ). This program starts with one of the 
unknowns found in step (i) above and sets its value equal to 1. Then it 
determines the maximum number of other unknowns that can be de- 
termined from this recursively, i.e., without having to invert any matrix. 
The best unknown to fix, the one that minimizes the number of un- 
knowns to be solved by inverting a matrix, is selected and set equal 
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to 1. The matrix corresponding to the equations for the other unknowns, 
and the right-hand sides for these equations, are then generated in 
symbolic form. The output of this program is then a subroutine which 
generates the coefficients of c r (s) in symbolic form, as functions of the 
given probabilities. Some excerpts are shown in Tables I and II. In Table 
I is the "triangular" part of the equations for k = 5. /2,- . . . ; 6 here is 
P, . . . ,- 5 normalized such that Poooooo = 1 and P(i) = Pr\x n = i - 1}. In 
Table II are equations for generating the coefficients of c r (s) from 

VI. COMPUTATION OF MARGINAL DISTRIBUTIONS 

Equations (37) and (38) give expressions for the generating functions 
of the equilibrium distributions of y>j n , j = 0,1, . . . ,k. As in (30), yon is 
b n the queue length. Hence, the equilibrium queue size distribution has 
as its generating function </> (s). If *oj = lim,,^ Pr{5„ = ;'}, ;' = 0,1, 
then 0o(s) = 2J=o ttojsJ. One way to find \ttoj\ is to start with #/j(s) and 
iterate using (38), thus obtaining an expression for (f>o(s), then to inverse 
transform o (s). For example, using s = e~i a , we can treat 0o(s) as a 
function of co, and then finding (7r ;} corresponds to finding the Fourier 
series for </>o(e~- ,w ). 

We will present a different method here. Generally, the quantities of 
interest are 

n , = £ tto, = lim Pr)6 n < ; | for ; < N, 



i = 



nt- 



Table I 



ftinill = ((1-P(1 ))IP( l))*#oooooo 
#111112 = (P(2)/P(l))*floooooo 

#222223 = (P(2)/PU))*#,Iim 
#222224 = (P(3)/P(2))*#,,111> 
#333335 = (P(3)/P(2))*ft 22 2223 
#333336 = (P(4)/P(3))*P. 222224 
#444447 = (P(4)/P(3))*fl 3 33335 
#444448 = (P(5)/P(4))*/?333336 
#000011 - P(l)*/?111112 
#111123 = <P(2)/P(1))*#000011 
#111133 - P(D*#222224 
#222235 = (P(3)/P(2))*#U1123 
#222245 = (P(2)/P(1))'# 11 1133 
# 222255 - P(D*#333336 
#333347 = (P(4)/P(3))*«222235 
#333357 = (P(3)/P(2))*P 2 22245 
#333367 =(P(2)/P(l))*#22225f. 
#333377 - P(l)*#444448 
#000122 = P(1)*#11U23 
#000222 = P(l)*#111133 
#111234 = (P(2)/P(l))*#000122 
#111334 = (P(2)/P(l))*#O0O222 
#111244 = P(D*#222235 
#111344 = P(D*#222245 
#111444 = P(l)*#222255 
#222346 = (P(3)/P(2))*# U ,234 



#222356 = (P(2)/P(1))*P.„1244 
#222446 = (P(3)/P(2))*#n,334 
#222456 =(P(2)/P(D)*#111344 
#222556 =(P(2)/P(D)*#111444 
#222366 — P(D*#333347 
#222466 = P(D*#333357 
#222566 = P(D*#333367 
#222666 = P(l)*#333377 
#001233 = P(l)*#111234 
#001333 = P(l)*#111244 
#002233 = P(l)*#111334 
#002333 = P(l)*#111344 
#003333 = P(l)*#111444 
#112345 = (P(2)/P(D)*#001233 
#112445 = (P(2)/P(l))*#001333 
#113345 = (P(2)/P(l))*#002233 
#113445 = (P(2)/P(l))*#002333 
#114445 = (P(2)/P(l))*#O03333 
#112355 - P(l)*#222346 
#112455 = P(l)*#222356 
#112555 = P(l)*#222366 
#113355 = P(D*#222446 
#113455 = P(l)*#222456 
#113555 = P(l)*#222466 
#114455 = P(l)*#222556 
#114555 = P(l)*#222566 
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Table II 



"000000 + #000011 + #000111 + #000122 + #000222 + #001111 + #001122 
C(l,l) + #001222 + #001233 + #001333 + #002222 + #002233 + #002333 + #003333 
#011111 + #011122 + #011222 + #011233 + #011333 + #012222 + #012233 
C(l,2) + #012333 + #012344 + #012444 + #013333 + #013344 + #013444 + #014444 
#022222 + #022233 + #022333 + #022344 + #022444 + #023333 + #023344 
C(l,3) + #023444 + #024444 
#033333 + #033344 + #033444 + #034444 
#044444 

#000000 + #000011 + #000111 + #000122 + #000222 

#001111 + #001122 + #001222 + #001233 + #001333 + #011111 + #011122 
C(2,2) + #011222 + #011233 + #011333 

#002222 + #002233 + #002333 + #012222 + #012233 + #012333 + #012344 
C(2,3) + #012444 + #022222 + #022233 + #022333 + #022344 + #022444 
#003333 + #013333 + #013344 + #013444 + #023333 + #023344 + #023444 
C(2,4) + #033333 + #033344 + #033444 
}«•»"/ - #014444 + #024444 + #034444 + #044444 

C(3,l) = #000000 + #000011 

C(3,2) = #000111 + #000122 + #0 

C(3,3) = #000222 + #001222 + #0 



C(l,l) - #000000 + # 

C(l,l) 

C(l,2) 
C(l,2) 
C(l,3) 
C(l,3) 
C(l,4) 

C(l,5) - #044444 

C(2,l) " 

C(2,2) 

C(2,2) 

C(2,3) 

C(2,3) 

C(2,4) 

C(2,4) 

C(2,5) 




^011233 

C(3,3) 

C(3,4) 
C(3,4) 
C(3,5) 
C(3,5) 



C(3,3) + #012222 + #012233 + #022222 + #022233 

#001333 + #002333 + #003333 + #011333 + #012333 + #012344 + #013333 

C(3,4) + #013344 + #022333 + #022344 + #023333 + #023344 + #033333 + #033344 

#012444 + #013444 + #014444 + #022444 + #023444 + #024444 + #033444 

C(3,5) + #034444 + #044444 

#(KMXKXI 



>-W»«/ — ^(u>"l ' «»U34444 i • MJ44444 

C(4,l) = #000000 

C(4,2) = #000011 + #000111 + #001111 + #011111 

C(4,3) = #000122 + #000222 + #001122 + #001222 + #002222 + #011122 + #01 

C(4,3) = C(4,3) + #012222 + #022222 

C(4,4) = #001233 + #001333 + #002233 + #002333 + #003333 + #011233 + #011333 
C(4,4) = C(4,4) + #012233 + #012333 + #013333 + #022233 + #022333 + #023333 + #l 
C(4,5) = #012344 + #012444 + #013344 + #013444 + #014444 + #022344 + #022444 

: C(4,5) + #023344 + #023444 + #024444 + #033344 + #033444 + #034444 + # 



1222 

'033333 



C(4,5) - #012344 + #012444 + #013344 + #013444 + #014444 + #022344 + #022444 

C(4,5) = C(4,5) + #023344 + #023444 + #024444 + #033344 + #033444 + #034444 + #044444 

where N is some initially selected constant. The method presented here 
determines IIo;, j < N, explicitly in a finite number of additions and 
multiplications. Let r (s) = 2"=o itrjsJ and H r j = 2/ =0 ir r i. We can write 
(37) as 

0fc(s)= : — 77 • (24) 

(pkvKS) ~S 

Denote </>*(s)/(l — s) by fais), so that iM«) = 2 f=o RkjS* for \s\ < 1. 
Therefore, 

Ms) = , |s|<l. (25) 

(Pkv(S) - S 

Similarly defining ^y(s) = 2° =0 n r /»- / for r = k — l,k — 2, . . . ,0, we 
have 

Ms) = *- x [*r+i(») " c r+ i(s)]0 ry (s). (26) 

The functions ry (s) are, of course, determined from the distributions 
of x„ and the constants a). 

For each r, the process of determining U. r j from Uh-ij can be described 
as follows. Subtract the known constants c r+ ij from n r +ij fory < degree 
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5rj = 



of c r+ i(s). Since n r+ io = c r +i,o, the sequence b r j corresponding to 
s -1 [i/v+i(s) — c r+ i(s)] is such that b r j = for,/ < and 

n r+ ij+i - Cr+ij+i for 1 < j + 1 < degree of c r+1 (s) 

n r+ i,y+i for; + 1 > degree of c r +i(s) 

In order to get {n r; ), we now convolve this sequence \5 r j\ with the se- 
quence corresponding to rL ,(s), say \p r j\- Therefore, n r; = S/=o 5 r j-,p r ,-. 
This process involves only a finite number of multiplications and addi- 
tions as long as only a finite number of II r /s are sought. Notice that in 
order to determine Iloy, j < N, we have to start with values of Rkj for j 
^N + k. 

We will now show that there is a recursion which allows us to compute 
Hjy, j < N + k, in a finite number of arithmetic operations. Let IT*, 
correspond to <f>kv(s)/(<t>kv(s) ~ s ) and let 0fc y (s) = 2" =0 PkjS*. Then, 
equating coefficients of like powers of s on both sides of 

>=0 <Pkv\S) - S 

we can derive the following. 

II'*o = 1, 

n'kj = (Pkj + n'v-i - i PkiWk,i-i)/Pk<>, j = 1,2, .... (28) 

i=l 

Once the U'kj have been determined, we convolve the sequence {11'*/} 
with [ckj\ to get {II kj], as seen from (25). Summarizing, we have shown 
that Iloy for; < N can be determined from c r (s), r=l,...,k and <f> ru (s), 
r — 0,1, ... jk y by performing only a finite number of multiplications and 
additions. The method described was used in calculating the probabil- 
ities presented in Section VII. 

VII. AN EXAMPLE OF THE CALCULATIONS 

We calculated the queue size distributions for various traffic inten- 
sities Ez n for the following queuing models: 

(i) b n+1 = (b„ - 1) + + 2x„, k = 

(ii) b n+ i = (b n - 1)+ + x n + x n - 3 , k = 3 

(Hi) b n +i = {b n - 1) + + x n + x n - 4 , k = 4 

{iu) b n +i = (b n - 1) + + x n + x n - 5 , k = 5, 

where the i.i.d. random variables x n are assumed to be distributed ac- 
cording to the Poisson law. These cases are referred to equivalently by 
referring to the value of k . The objective was to determine what buffer 
size would suffice for a concentrator used to buffer terminals that gen- 
erate two packets per message but are slower than the trunk line. We 
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present here some of the results. In Fig. 3 the abscissa corresponds to 
traffic intensity: average number of packets transmitted on the trunk 
line per unit of time. In the steady state, the probability that the queue 
size exceeds that shown on the ordinate is less than 10 -4 for each value 
of traffic intensity. The cases when all packets arrive at once {k = 0) or 
"infinitely apart" (k \ °°) are shown, as well as the case k = 5. When k 
= 0, the queue size corresponding to each value of traffic intensity is 
either twice that of the case k = « , or one less than twice. 4 These two 
cases are further compared in Fig. 4, this time for buffer sizes 20 and 40, 
and the logarithm to base 10 of the probability of the queue size ex- 
ceeding 20 and 40 is plotted. 

For fixed traffic intensities, the change in the probability of overflow 
as a function of buffer size can be seen from Figs. 5, 6, and 7. From these 
figures we can see that for low traffic intensities and low buffer sizes there 
is a difference between the batch case k = 0, and k = 5 (see Fig. 5), but 
for larger traffic intensities the difference decreases substantially. From 
the formulas (37) and (38), it can be shown that the tail of queue size 
distribution is geometric for each value of k. Furthermore, for all finite 
values of k, the common ratio is the same as that of the case when k = 
0, so that the similarity in the behavior of queue size distributions for 
large values of queue size is as expected. The slopes of the curves marked 
k = 5 in Figs. 5, 6, and 7 approach those of the curves marked k = for 
large values of the abscissa. Hence, in applications where probabilities 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.E 

TRAFFIC INTENSITY 



0.9 1.0. 



Fig. 3 — Buffer size vs traffic intensity for k = 0, 5 and <=, and probability of overflow 
< 10" 4 . 
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Fig. 4 — Probability of overflow vs traffic intensity for k = and °°, and buffer sizes 20 
and 40. 

of overflow of 10~ 7 or smaller are required, the relative slowness of the 
sources of packets does not seem to reduce the buffer size required. 
Packets may be assumed to arrive simultaneously for the purpose of 
estimating the buffer size. 



APPENDIX A 
Summary of Formulas 

We here summarize formulas for calculating the equilibrium queue 
size distribution. With a model for the input process z n as in (2), the 
queue size b n is described by 



/ k 



b n+1 = (b n -l)++ E L «K-;. 

f=l./=0 



(29) 



Various formulas pertaining to (29) were derived. 4 The reader is referred 
to Ref. 4 for proofs of the formulas presented here. Define yo n = b n and, 
for r = 0,1, . . . ,k - 1, 



ik 

y r +i,n = y™ + L L awn-j+r- 

i'=l ;=r+l 



(30) 
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5 10 

BUFFER SIZE 



Fig. 5 — Probability of overflow vs buffer size for k = 0, 5 and », and traffic intensity 
0.4. 

Note thatyon <yin < • • • < ykn- The vector process (yo n .yin, • • • dkn), 
n = 0,1,2, . . . , is Markov under the assumption that (x£,*i, ■ ■ • ,x l n ), n 
= 0,1,2, . . . , is a vector sequence of independent identically distributed 
random variables. The Markov chain, denoted by S, which corresponds 
to (yon.yin, ■ ■ ■ ,ykn), was shown to be positive recurrent when Ez n < 
1. The states of S are those which communicate with (0,0, ... ,0), since 
it is assumed that the buffer is empty and no one is transmitting initially, 
at n = 0. Let 



L «} = Mr, L l»\X l n = V m , r = 0,1, . . . ,k. 
7=0 1=1 



(31) 



The probabilities that enter the calculations turn out to be only those 
corresponding to the random variables u rn , r = 0,1, ... ,k. Let pi 0i i h . . . j k 
= Pr{t; n = io» u in = h, • ■ ■ ,Vkn m ik\- Then the transition probabilities 
for S are given by: 



pn +1 = Y n . ■ . . P? 

J0=0,jtA 



+ E pMTVi+l.ii-72+l, 

jo>0.UA 



,ik-jk+iPy ( 32 ) 
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10 20 

BUFFER SIZE 



Fig. 6 — Probability of overflow vs buffer size for k = 0, 5 and °° , and traffic intensity 
0.6. 

The sums in (32) are over those indices JoJi, . . . Jk which correspond 
to states communicating with (0,0, ... ,0) denoted by A. As mentioned 
above, when Ez n < 1, then lim„ta> P" = P\ exists and Pi satisfies (32) with 
Pi and P" +1 both replaced by Pi. The generating function corresponding 
to P h 



0(80,5!, . . . ,Sfe) - EPtfftep...*? 

UA 



satisfies 



(f)(s ,Si, . . . ,S k ) = [0(l,So,«b ■ • • ^k-^Sk-lSk) f[ s i 1 

i-0 

+ ( 1 - n s i _1 ) 0(O,S O ,Si, . . . ,S k -2,Sk-lSk)]<t>v(S0,Si, . . . ,s k ), 

where 

k 
faisoM, ...,s k ) = E n sf in . 

i=0 

It can be shown 4 that <j> has the representation 

<f>(s ,si, . . . ,s k ) = £ Cjdj(s , . . . ,Sk), 
i 



(33) 



(34) 
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TRAFFIC INTENSITY = 0.8 




6 



K = ~ 



10 20 

BUFFER SIZE 



30 



Fig. 7 — Probability of overflow vs buffer size for k = 0, 5 and °°, and traffic intensity 
0.8. 

where 0j(sn, . . . ,Sk) are solutions of ' 

k 

Oj(so,Si, . . . ,S k ) = [0j(l,S O ,Sl, • • • ,Sk-2,Sk-lSk) II Si -1 

+ M (l - I! sr^W 2 • • • si k --l(s k -iSk) jk ](t>o(so,si, . . . ,s k ), (35) 

j=0 

with n = 1 — Ez n = 1 — Evkn- Substitution of (34) into (33) yields a sys- 
tem of equations for c } which together with 2cj = 1 uniquely deter- 
mine 

0(O,si, . . . ,8 k ) =Ml c iS { 1 . . . s{ k . (36) 

j 

Once 0(O,si, . . . ,Sk) is found, the generating functions corresponding 

to the marginal distributions, namely lim n f «, Es yin = 0,-(s), are solutions 

of 



0Jfe(s) = 



(1 -S ^CkJS^kuis) 
1 -S -1 0*u(s) 



(37) 



where 0/y(s) = Es w > B , J = 0,1, . . . ,k, and, for r ■ k — 1, . . . ,1,0, 

<Ms) = [s-Vr+its) + (1 - s- 1 )c,+i(s)]0 ri ,(s). (38) 
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Here the polynomials c r (s), r = 1, . . . ,k are of finite degree and are given 
by 

r 

c r {s) = 0(0,1, . • • ,s,l, • • • ,D = M L C,L,2 . . . j r . . . Ik* lr - ( 39 ) 

i 

Then (37) uniquely determines <t>k(s). Using (38) k times yields <f>o(s), 
which by definition of y 0n [see (30)] is the generating function of 
steady-state queue size. 

An alternate representation of the generating function corresponding 
to Pi is obtained by setting Uj = 11*=, s it j = 0, . . . ,k and defining 

0(s o ,si, • • • ,Sk) = *("o."i, • • • Mk)- (40) 

Then, corresponding to (33), 
$(t*o,Ui, . . . Mk) = [wo 1 ^(wo»"o»wi. • ■ • Mk-i) 

+ (1 - Uo 1 )*(0,U ,Wl, • • • ,U fc -i)]* w (uo,Ui, . . . ,Ufc), (41) 

where 

*„(U ,M1, • • • ,"*) - E ( ft <'"), "rn = L «{*'„. (42) 

It follows from (41) that, for; = 0,1, . . . Jk - 2, (k > 2), 

$(s, . . . ,s,ui, ■ ■ ■ Mk-j) = [s _1 $>(s, . . . ,s,u\, ■ • ■ Mh-j-l) 

+ (i - 8 _1 )*(0a . . . ,s,m, ■ • ■ Mk-j-i)] 

X $„(s, . . . ,s,ui, ■ ■ ■ Mk-j), (43) 
and 

$(s ( . . . ,s,m) = [s" 1( l>(s, . . . ,s) 

+ (1 - s- x )*(0,s, . . . ,s)]4> u (s, . . . ,s,ui). (44) 

If we set ui = s in (44), and solve for <i>(s, . . . ,s), we obtain 

. (1 -g)$(0,s,. . . ,s)* u (s, . . . ,s) . 

[*„(«, . . . ,s) -s\ 

It was shown 4 that $(0,Ui, . . . ,u*) is a multinomial independent of 
u fe . From (43) to (45), $(s,Ui, ■ ■ ■ Mk) may be expressed in terms of 
*(0,s, . . . ,s,Uh • • • Mk-j-i), ) = 0, . . . ,k - 2, and «i>(0,s, . . . ,s). If we let 
s — in this expression, and equate $(0,ui, .-.,"*) with the finite part, 
we obtain a system of homogeneous linear equations for the coefficients 
in the multinomial. In general, we also obtain a (consistent) set of ho- 
mogeneous linear equations from finiteness conditions. The normali- 
zation condition is *(0,1, • • ■ ,D = M- From (36) and (40), it follows 
that 

*(0,m, . . . Mk) = m Z c h . . . Jh u{ 1 n m^"*- 1 . (46) 

j i-2 
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Also, from (39), 

r k—r 

C r (f)-*(ta...Al,....D. (47) 

Formulas for calculating the first two moments of the equilibrium 
queue size distribution are derived in Appendix B. 

APPENDIX B 

First and Second Moments 

We here derive expressions for the first two moments of y r = lim n — «> 
y rn . We note that yo represents the equilibrium queue size. By defini- 
tion, 

r ( S ) = E( S yr), ry ( S ) = E( S "rn) ) ( 48 ) 

for r = 0, . . . ,k . Hence, 

r (l) = l, <f>' r (l) = Ey r , <t>" r (l)=E(y?)-Ey r (49) 

and 

6,(1) -1, <t>' n (l) = Ev m , (i>"r u a)=E(u? n )-Ev rn . (50) 
We also note that, from (39), since 2 c } = 1, 

c r (l) = n = 1 - Ev kn , r=l,...,k. (51) 

From (37), 

4ft<9)-C*(9)A*(«)/(s), fb) = . (1 " S) i ' < 52 > 

and, from (38), 

Ms) = [*"Vr+l(«) + (1 - S-^Cr+lC*)]^^), r = 0, . . . ,/z - 1. 

(53) 
If we differentiate (53), we obtain 

0',(S) = [-S- 2 r+1 (s) + S~ Vr+l(s) + S" 2 C r+ i(s) + (1 - S^c'^-lOO] 

x ru (s) + [s-^r+iOO + (1 - a-^br+iCaH^U). (54) 
If we set s = 1 in (54), and use (49) to (51), we obtain 

Ey r = Ey r+l + Ev rn - Ev kn , r = 0, . . . ,k - 1. (55) 

Next, if we differentiate (52), we find that 

0'*(«) = [<?'*(«)**„(*) + C*(«)0'fa,(*)]/(«) + c*(s)<^(s)/'(s). (56) » 
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But 

<t>kv(s) = 1 + (S - l)*'to(l) + - (S - l)Vto(«) 

+ 7<s- Wfa,(D+ .... (57) 
b 

Hence, from (52), 

,, x 1 , (s - DVkM 

f( s )=T, T-77TT + 



[1-0'Aud)] 2[l-0' fcu (l)] 2 

+ . . . . (58) 



, («-i) 2 ( £fafl) | fo'tod)] 2 



2 Ull-^tod)] 2 2[l-tf'fa,(D] 3 
If we let s -* 1 in (56), and use (49) to (51) and (58), we obtain 

E y k = Tx — p — t + Eukn + ~^\ — v — \ ■ (59) 

(1-EUkn) 2(1 -EVkn) 

From (55), it follows that 

E yj = Ey k + *£ Ev m - (k - j)Eu kn , j - 0, . . . ,k - 1. (60) 

This determines Eyy, and in particular Ey , in view of (59). 

For the second-order moments, if we differentiate (54) and set s = 1, 
we obtain 
E(y r 2 ) = E(yhi) ~ Ey r+ i + Ey r + E(v? n ) - Ev m 

+ 2(1 - Eu rn )(Ev kn - Ey r+1 ) + 2c' r+ i(l), (61) 
for r = 0, . . . ,k — 1. Finally, if we differentiate (56) and let s — I, and 
use the relationship 

0" to (D = E(vi n ) - mvL) + »y, (62) 

which follows from (48), we obtain, after some simplification, 
-, 2x „ , \c" k (l) + 2{Ev kn )c' k (l)} , [E^D-^fenKfed) 

[£(t;L)--g^n] [g(^D--g^] 2 ( g 3) 

Zd-EVkn) 2(1- EVkn) 2 

An expression for E(yo) may be obtained from k applications of (61), with 
the help of (59), (60), and (63). 

APPENDIX C 
Determination of Constants 

We first show how to determine which constants Cj occur in (46) for 
the example of (3). From (2) and (30), 

yOn = b n , y r +l, n - yrn = *n-k+r, T = 0, . . . ,k - 1. (64) 

Hence, 

"0" 11 U? l n-k+r-\). (65) 

r=l / 
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But from iteration of (29), it is evident that 

b n » -► £ Zn-j^ I - 1, l=h... ,k, (66) 

where z n is given by (2). Hence, for the example of (3), 

b n = => L *i-i < I - 1, J - 1, . . . ,*. (67) 

The inequalities in (67) determine the admissible vectors (x n - k , . . . ,x n -i) 
corresponding to b n = 0, and thus, from (65), which constants Cj occur 
in (46). Note, in particular, that *J_i = implies that $(0,ui, . . . ,Uk) 
is independent of Uk. 

For purposes of illustration, we show how to calculate the values of 
Cj for the example of (3), subject to (4), in the case k = 3. From the above 
procedure it is found that 

$(0,Ui,U2,U3) = X("l,"2) = M(C(K)0 + Ciulli + Cq U U 2 

+ C222"l + Cl22"l"2)- (68) 

But from (2) and (42), 

*i/("o,wi,"2,"3) = 6[(1 - p)u u 3 + pw ]^(w ). (69) 

Hence, from (45), 

{e[(i -p)s 2 + ps]^(s)-s} 

Then, from (44), we obtain 

,, v (1 - s) x (s,s)e[(l - p)su 1 + ps]*(s) mx 

$(s,s,s,ui) = — — — — 777- : • (71) 

{6[(1 - p)s 2 + ps]V(s) - 9) 

Also, from (43), we have 

$(s,s,Ui,u 2 ) = [s _1 $(s,s,s,wi) + (1 -s _1 )x(s,s)] 

X 9[(1 - p)su 2 + ps] ¥(a), (72) 

and 

*(s,ui,u 2 ,w 3 ) = [s _1 $(s,s,ui,u 2 ) + (1 -s _1 )x(s,"i)] 

X 9[(1 - p)su 3 + ps]Ms). (73) 

From (71) to (73), it follows that 

(1 — s) 

SfotfbU&Us) = ©[(1 - P)«"3 + ps]*(s) 

s 



X 



e[(i - p)su 2 + ps] — 

s 



r e[(i- P )su 1 + ps)ns) I 1 

X ^77, TTT — TT7-; : - 1 x(m) ~ x(«,"i) • (74) 

L|0l(l - p)s 2 + ps]V(s) - s\ J 
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From finiteness at s = we deduce, with the help of (6), that 

x(0,ui) = m[1 + (1 - p)P i<7o"i]cooo- (75) 

Hence, from (68), 

con = (1 - p)Pi<7oCooo- (76) 

If we now let s — in (74), and equate coefficients, we obtain, in addition 
to (76), the relations 

C122 = (1 - P) 2 P1 2 <?0C000» C222 = (1 ~ P) 2 PoP 2<7 0^000, (77) 

and 

cm = (1 - p)<7o(Pi + PoPiQi + 2ppoP2Qo)cooo 

+ (1 - p)poPi<7o(cm + con) - Po<7oCi22- (78) 

If we substitute the values of Con and C122, from (76) and (77), into (78), 
we obtain (9). 
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